Risk stratification of new-onset psychiatric disorders using clinically distinct traumatic brain injury phenotypes

Background Patients with traumatic brain injury (TBI) constitute a highly heterogeneous population, with varying risks for New-onset Psychiatric Disorders (NPDs). The objectives of this study were to identify TBI phenotypes and determine how NPDs differ among these phenotypes. Methods Hospitalized TBI patients from 2003 to 2019 were obtained from the provincial trauma registry. Propensity score matching was conducted to balance covariates among patients with TBI and controls. To uncover heterogeneity in TBI, latent class analysis (LCA)-based clustering was applied. LCA was conducted separately for two TBI cohorts: those with and without pre-injury psychiatric conditions The effect of classes on NPDs was assessed using log binomial regression models. Results A total of 3,453 patients with TBI and 13,112 controls were included in the analysis. In a conditional regression involving propensity matched patients with TBI and controls, TBI was significantly associated with the development of NPD-A (OR: 2.78; 95% CI: 2.49–3.09), as well as NPD-P (OR: 2.36; 95% CI: 2.07–2.70). Eight distinct latent classes were identified which differed in the incidence of NPDs. Four classes displayed a 53% (RR:1.53; 95% CI: 1.31–1.78), 48% (RR:1.48; 95% CI: 1.26–1.74), 28% (RR:1.28; 95% CI: 1.08–1.54), and 20% (RR: 1.20, 95%CI: 1.03–1.39), increased NPD risk. Conclusion TBI is a significant predictor of NPDs. There are clinically distinguishable phenotypes with different patterns of NPD risk among patients with TBI. Identifying individuals with respect to their phenotype may improve risk stratification of patients with TBI and promote early intervention for psychiatric care in this vulnerable population. Supplementary Information The online version contains supplementary material available at 10.1186/s13690-024-01346-w.


Introduction
Traumatic brain injury (TBI) is a significant contributor to global injury burden, with an age-standardized incidence rate of 369 per 100,000 population [1].Each year in Canada, there are approximately 166,455 TBIs [2], resulting in roughly 20,000 hospitalizations and annual direct medical costs of $120.7 million CAD [3].TBI is a leading cause of mortality in Canada, accounting for approximately 23% of all injury-related deaths [2].
New-onset psychiatric disorders (NPDs) are common after TBI [4].Factors such as the location of brain lesion, prior lifetime psychiatric disorders, family function, family psychiatric history, and socioeconomic status have shown significant associations with the development of NPDs following pediatric TBI [5][6][7][8][9].Furthermore, NPDs are more common in adults with a history of pediatric TBI compared with controls without childhood head injuries, and the presence pre-injury lifetime psychiatric disorders notably predicts the occurrence of NPDs [10].
Despite the growing body of research on the link between TBI and NPDs [4], there remain significant gaps in our knowledge.Pre-injury psychiatric history is a wellestablished risk factor for NPDs in patients with TBI [7,10,11].Several studies have explored the risk of NPDs in patients without a prior psychiatric history, finding that many such patients also develop post-injury disorders [12,13].Yet few studies have compared the risk of NPDs in TBI patients with and without pre-injury psychiatric disorders.Secondly, TBI patients constitute a highly heterogeneous population, with different profiles of risk and long-term outcomes [14][15][16].While most studies have examined the relationship between various predictors and NPDs following TBI, none have explored the heterogeneity in the relationship between demographics, injury, psychiatric history and NPDs in patients with TBI.
The current study addresses these knowledge gaps by defining NPD as two distinct outcomes experienced by those with and without a psychiatric history.A large population-based dataset of patients with TBI and controls was used to demonstrate the difference in the development of NPDs over two years.Using demographic, injury, and clinical history data, clustering models were specified and fitted to identify phenotypes of patients at high risk of NPDs.We hypothesized that this approach could identify patient clusters, each with a distinct profile of clinical characteristics associated with significant differences in NPDs.
This investigation addressed three specific research questions: (1) What is the effect of TBI on the development of NPDs in a propensity-matched sample of injured patients and population-based controls?(2) Within the TBI cohort, are there specific phenotypes which differ with respect to demographics, injury variables, and preinjury psychiatric conditions?(3) Are these phenotypes associated with the development of NPDs?

Data sources Nova Scotia Trauma Registry (NSTR)
The TBI cohort was obtained from the NSTR, a comprehensive, provincial, population-based trauma registry.NSTR is an extensive database that captures demographic and clinical data on all major trauma cases requiring hospitalization in the province.

Canadian Institute of Health Information-Discharge Abstract Database (CIHI DAD)
CIHI DAD captures administrative, clinical and demographic information on hospital discharges (including deaths, sign-outs and transfers).Psychiatric diagnoses were based on the International Classification of Diseases, Tenth Revision (ICD-10) codes F04-F99.

MSI Physician Billing (MED)
MED records physician billing information with details of a service encounter between an individual and a provider.Psychiatric diagnoses were based on the International Classification of Diseases, Ninth Revision (ICD-9) codes 210-319.

National Ambulatory Care Reporting System (NACRS)
NACRS contains data for all hospital-based and community-based ambulatory care.Psychiatric diagnoses were based on ICD-10 codes F04-F99.

Exposure
The exposure of interest in this study was the occurrence of TBI.We included patients with TBI admitted to a provincial hospital between 2003 and 2019.The severity of TBI was determined using the Abbreviated Injury Scale (AIS).Patients presenting solely with scalp abrasions or superficial lacerations (not representative of intracranial injury), are often categorized as AIS Head score = 1 [17], and were excluded from the study cohort.Only subjects with a maximum AIS Head score ≥ 2 were included in the study.This definition for TBI severity includes a broad range of injury severity and has been used in previous studies [17][18][19][20].Controls were obtained from the provincial health insured registry and did not have a history of TBI or develop TBI in the two years of follow up based on clinical encounters (See Supplementary Data).

Covariates
Guided by the literature, age, sex, medical comorbidities (defined through the Charlson Comorbidity Index [CCI]), pre-injury psychiatric diagnoses, AIS Head score, injury mechanism, discharge destination and length of stay (LOS) in acute care facility were identified as important covariates to consider as potential confounders.

Outcome
The outcome of interest was NPD in the first two years post-TBI.Previous studies have defined NPD (also referred to as novel psychiatric disorder) to occur in one of two conditions [6][7][8]: (1) the development of a psychiatric disorder after injury in a patient with no lifetime preinjury psychiatric disorder; or (2) the development of a post-injury psychiatric disorder in a patient with a lifetime psychiatric disorder that was never before present (e.g., a subject with a lifetime history of major depressive disorder who develops anxiety disorder after the injury would receive the classification, but would not if only a new episode of major depression occurred).In many cases, the "lifetime psychiatric disorder" variable refers to any psychiatric disorder present prior to the injury [10,13].In the current study, we defined NPD as two distinct outcomes, namely, NPD in the absence of psychiatric history (NPD-A) and NPD in the presence of psychiatric history (NPD-P).These outcomes were defined to understand influence of pre-injury psychiatric history on the occurrence of NPDs.

Study design
The overall study design is presented in Fig. 1.The study is a matched cohort design with a two-year followup for outcome collection.The study unfolded in two phases: firstly, we assessed the development of NPDs in a propensity-matched cohort of patients with TBI and population-based controls; secondly latent class analysis (LCA) was applied to the TBI cohort, identifying TBI phenotypes.

Propensity score matching
The covariates entered into the propensity score were age, sex, CCI group (no comorbidities [0], mild [1][2], moderate [3][4], and severe [≥ 5]), and psychiatric diagnoses (ICD9:290-319, ICD10: F04-F99).Cases and controls were matched using greedy nearest neighbor matching [21] without replacement using a calliper width of 0.2 standard deviations of the logit of the propensity score.A matching ratio of 4-to-1 (up to 4 matches of controls per case) was used (See Supplementary Data).The matched cohort was then divided into those with and without preinjury psychiatric conditions before further analysis.

Latent class analysis
LCA was used to derive clinical phenotypes of patients with TBI, which represent their distinct underlying individual attributes.Details on the LCA methodology are described in the Supplementary methods.Since preinjury psychiatric disorders are common among individuals with TBI, we aimed to study the influence of TBI on those with existing psychiatric disorders as well as those without pre-injury psychiatric conditions; LCA was conducted separately for these two cohorts.Confirmation of ICD codes related to psychiatric conditions in the two years prior to the injury were used to define pre-injury psychiatric conditions.For those without psychiatric history, indicators included sex, categories of age, categories of CCI, injury mechanism, injury severity (maximum AIS Head), prolonged LOS, and discharge destination (Table 1).For those with pre-injury psychiatric disorders, we used the same indicators and added past psychiatric conditions as additional indicators.Only psychiatric conditions with a prevalence of > 10% were included: disorders relating to sleep, organic conditions, mood, anxiety, drug substance abuse disorder (SUD), and alcohol SUD.For age, we utilized four categories representing different stages of life: pediatric (0-18 years), young adult (19-35 years), middle-aged adult (36-59), and older adult (≥ 60 years).LOS was dichotomized as brief or prolonged, defined as ≥ 21 days based on previous research [22,23].Relative fit was evaluated using the Bayesian Information Criterion (BIC), the sample-size adjusted BIC (SABIC), and Vuong-Lo-Mendell-Rubin adjusted likelihood ratio test (VLMR-LRT) [24].Empirical robustness was complemented by interpretability during the process of model selection [25].
Characteristics were compared across the identified phenotypes.Categorical variables were reported as number (percentage) and were compared across classes using the chi square test.In the overall matched cohort, conditional regression was used to account for the matched sets and to provide odds ratios (OR) as the measure of association between the exposure (TBI) and the outcome

HDNS 2
Pre-injury psychiatric conditions 3   Sleep disorders, Organic disorders, mood disorders, anxiety disorders, drug-related substance abuse, alcohol-related substance abuse.Comorbid psychiatric conditions were defined as > 1 pre-injury psychiatric conditions.(NPD) [26].Following the identification of latent classes, a conditional logistic regression analysis was also conducted within each phenotype to compare the odds of NPD in TBI patients with their matched controls.The OR from conditional regression estimates the relative odds of NPD occurring in the TBI group compared to the control group, while controlling for the matching variables.Log-binomial regression models were fitted to observe the association between phenotypes and the outcome of NPD.We used log-binomial regression to directly estimate the relative risk (RR), which is more interpretable than the OR when the outcome prevalence is high (commonly defined as > 10%) [27].A p < 0.05 was considered as statistical significance.

Characteristics of study population
There were 3,453 patients with TBI and 13,112 matched controls in the analysis.Table 2 shows the frequencies for the matched variables between cases and controls.TBI and control groups did not significantly differ in age, sex, comorbidities, or history of psychiatric disorders.In the two-year follow-up, a significantly greater proportion of TBI patients developed NPD-A (35.6% vs. 16.9%,p < 0.001) and NPD-P (46.9% vs. 27.7%,p < 0.001) than matched controls.Psychiatric diagnostic categories are described in Supplementary Table 1.

Comparison of TBI vs. non-TBI patients for NPDs
In conditional regression analysis performed on the matched groups in the overall study population (Table 3), TBI was significantly associated with the development of NPD-A (OR: 2.78; 95% CI: 2.49-3.09),as well as NPD-P (OR: 2.36; 95% CI: 2.07-2.70).

Latent class analysis
Table 4 presents fit indices for latent class models ranging from one to five classes.A 4-class solution provided the best fit for both TBI cohorts.Supplementary Tables 2-3 provide the mean posterior probabilities of the selected models, representing the average likelihood of individuals being assigned to their respective latent classes.All class assignments were > 80%, demonstrating high class assignment certainty for individuals.

Class descriptions for TBI cohort without pre-injury psychiatric conditions
Table 5 summarizes the descriptive statistics of the four classes and a visual representation is shown in Fig. 2. Significant differences were observed across all characteristics between classes.Additionally, the incidence of NPD-A showed a significant variation between classes (p < 0.001).

Older adults/extended recovery/supportive care class
This class was composed largely of older adults (98%) who experienced falls (89%).They were discharged to other acute care facilities (43%) or home with support services (29%) after a prolonged LOS (46%).

Young males/extended recovery/rehabilitative care class
Members of this class were young adult (42%) and middle age (38%) males (82%).Their head injury was mostly caused by MVCs (69%) and resulted in a severe injury (AIS = 5; 42%).The majority of this class had a high probability of being discharged to a rehabilitation facility (70%) after a prolonged LOS (84%).Within each of the identified phenotypes, the odds of NPD-A were significantly higher among TBI patients than their matched controls.The risk ratio of NPD-A in the older adult/extended recovery/supportive care class was 53% higher than the young males/timely discharge/ home recovery class (RR:1.53;95% CI: 1.31-1.78,Table 5).Additionally, the young males/extended recovery/ rehabilitative care class exhibited a 48% higher RR (RR:1.48;95% CI: 1.26-1.74)and the older males/timely discharge/home recovery class demonstrated a 20% higher RR (RR: 1.20, 95%CI: 1.03-1.39)than the young males/timely discharge/home recovery class.

Class descriptions for TBI cohort with pre-injury psychiatric conditions
Table 6 summarizes the estimated profiles of the 4-class model for the TBI cohort with pre-injury psychiatric conditions and a visual representation is presented in Fig. 3. Significant differences were observed across all indicators between classes.Moreover, the incidence of NPD-P showed a significant variation between classes (p = 0.02).The distinct classes identified in this cohort are as follows:

Older females/intermediate psychiatric burden class
Members of this class were older (98%) females (60%) with a high probability of falls as the cause of TBI (94%).This class was discharged home with support services (23%) or to a nursing home (30%).These individuals had pre-injury organic disorders (59%), mood disorders (29%), anxiety (27%), and a 49% probability of comorbid conditions.

Anxiety predominant/low psychiatric burden class
This class was composed of middle-age (46%) and older adult (53%) males (69%) who experienced falls (71%).They were discharged home (64%).Pre-injury anxiety (46%) was the most common disorder, but there were but no comorbid disorders in this class.
Within each of the identified phenotypes, the odds of NPD-P were significantly higher among TBI patients than their matched controls.The psychiatric complexity/ Suppression of cell counts: If any cell is < 5, the value (and its corresponding percentage) is suppressed and indicated with dash (-).If a cell is < 5 and only one value is suppressed in a row or column, the next highest value in that row or column is also suppressed 2 NPD-A was not used as an indicator in LCA.Differences in NPD-A incidence were determined after classes were modeled REF: reference level *p-value < 0.05 high comorbidity class exhibited a 28% higher risk (RR:1.28;95% CI: 1.08-1.54) of NPD-P compared to the young adult/low psychiatric burden class (Table 6).

Discussion
In the present study, a matched cohort was used to determine the increased risk of NPDs in those with and without TBI.Additionally, the potential effects of demographics, injury variables, medical comorbidities, and pre-injury psychiatric conditions on the NPDs experienced by TBI patients were described.We established that TBI is a heterogenous population demonstrating distinct clinical profiles (phenotypes), that show significant differences in the development of NPDs, suggesting their potential utility in risk stratification.This study is a proof-of-concept and further research is needed to assess the utility of these phenotypes to guide clinical management of TBI patients.
Previous research indicates that patients with TBI are at risk for developing NPDs with rates ranging from  Table 6 Analysis of patient characteristics by latent class for the TBI cohort with pre-injury psychiatric conditions.Adjusted relative risk (RR) of NPD-P is estimated from log binomial regression model.Odds ratio (OR) comparing patients with TBI to their matched controls is estimated using conditional regression for matched subjects 18.3 to 60.8% in the first year post-injury [11,28].Our results are consistent with these findings as the rates of NPD-A and NPD-P at two years in the current study were between 36 and 47%.Results from a conditional regression model in the matched cohort revealed that in the cohort without pre-injury psychiatric conditions, the odds of developing NPD were 2.8 times greater in those with TBI.In the cohort with pre-injury psychiatric conditions, those with TBI had 2.4 greater odds for developing NPD compared to controls.Thus, head injury alone and in combination with pre-injury psychiatric conditions appears to contribute to increased risk of post-injury psychiatric disorders compared to controls.A previous population-based study demonstrated an increase in risk for all psychiatric outcomes after head injury in patients without pre-injury psychiatric history [29].Prior studies also support the finding that among TBI patients, preinjury disorders are significant predictors of post-injury disorders [12,30].We speculate that there are several justifications for how TBI may act independently and in conjunction with pre-injury psychiatric disorders to increase the risk of NPDs.These mechanisms involve an interplay of disrupted connectivity, neuroinflammation, neurochemical imbalances, psychosocial factors, and environmental influences.Firstly, the neuroinflammatory response triggered by TBI can persist long after the initial injury.Chronic neuroinflammation has been linked to alterations in brain function and structure, affecting regions critical for emotional regulation, such as the prefrontal cortex and the limbic system, thereby contributing to the development and exacerbation of psychiatric disorders [31].Additionally, TBI affects the brain's structural and functional connectivity, leading to impairments in the communication between different brain regions involved in mood regulation, cognition, and behavior [32].Disrupted connectivity in the default mode network, salience network, and fronto-limbic circuits has been associated with various psychiatric disorders, including depression, bipolar disorder, schizophrenia, autism, and PTSD [33].TBI can lead to imbalances in neurotransmitters such as serotonin, dopamine, and glutamate, which play crucial roles in mood regulation, and their disruption can contribute to the development of various psychiatric disorders.Furthermore, the psychological distress of sustaining a TBI, such as loss of independence, changes in social roles, and difficulties with daily activities, may contribute to the development of psychiatric disorders like depression, anxiety, and PTSD [34].Finally, pre-existing genetic vulnerabilities and environmental factors, including stress, trauma, and lack of social support, may interact with the neurobiological changes caused by TBI, increasing the risk of manifestation of psychiatric conditions [4].
Among TBI patients without any pre-injury psychiatric history, the emerging phenotype was predominately determined by LOS and discharge destination while accounting for the effects of age-related vulnerability to psychiatric disorders.The young males/timely discharge/ home recovery phenotype was characterized by a notable proportion of young adults and pediatric patients and showed the lowest NPD-A incidence among all phenotypes.Age-related resilience against physical trauma and NPDs likely underpins the lower incidence of NPDs in this particular phenotype.The highest risk of NPD-A was evident in the older adults/extended recovery/supportive care class.Previous research has shown that prolonged LOS, especially in the intensive care unit, is a risk factor for long-term psychiatric disorders.Recovering at home with support services and relocation to a nursing home is a challenging transition for older adults, often resulting in high rates of persistent depression and anxiety [35].These environmental changes and age-related vulnerability coupled with the recovery from a head injury may exacerbate mental health symptoms with subsequent development of NPD-A.The young-middle-aged males/extended recovery/rehabilitative care phenotype was also associated with a high incidence of NPD-A.Suppression of cell counts: If any cell is < 5, the value (and its corresponding percentage) is suppressed and indicated with dash (-).If a cell is < 5 and only one value is suppressed in a row or column, the next highest value in that row or column is also suppressed 2 NPD-P was not used as an indicator in LCA.Differences in NPD-P incidence were determined after classes were modeled REF: reference level *p-value < 0.05

Table 6 (continued)
Interestingly, the severity of head injuries observed within this class was especially higher when compared to the other subgroups, suggesting that the extent of functional loss resulting from these injuries necessitated specialized rehabilitative care.Importantly, non-routine discharge to a supportive facility often implies that an individual has experienced functional and/or cognitive decline, which may serve as mediators in the emergence of NPDs.Individuals experiencing limitations in their ability to perform activities of daily living may develop feelings of helplessness and frustration, potentially leading to the development of NPDs.Similarly, cognitive decline after TBI may impede one's ability to process and cope with stressors, increasing vulnerability to psychiatric disorders.Therefore, non-routine discharge to supportive facilities serves as an indicator of functional and/or cognitive loss, which in turn can play a mediating role in the development of NPDs.The current study and data sources did not allow for an analysis of how functional, cognitive, and social outcomes of head injury may Fig. 3 Radar plots for pre-injury psychiatric burden for four phenotypes of TBI patients with psychiatric conditions be confounders, mediators, and/or modifiers of the association between TBI and NPDs.We acknowledge that a mediation analysis to quantify the extent to which the relationship between TBI and NPDs can be explained by one or more intermediate variables is critical in expanding our understanding of the complex causal pathway between head injury and NPDs.
We identified four distinct latent classes among TBI patients with pre-injury psychiatric history, each exhibiting a unique profile and highlighting the interaction between psychiatric burden and NPD-P incidence.The young adult/low psychiatric burden phenotype was characterized by the lowest burden of pre-existing psychiatric conditions and the lowest NPD-P incidence among all phenotypes.The reduced burden of psychiatric conditions is a likely explanation of the lower incidence of NPD-P in this particular phenotype.In contrast, the psychiatric complexity/high comorbidity class had the highest pre-injury psychiatric comorbidity levels and were discharged home post-injury.The high burden of psychiatric conditions coupled with the potential lack of screening and monitoring for NPDs in a home environment, may explain why this phenotype had the greatest NPD-P incidence among four classes.
Results of conditional logistic regression for each phenotype allowed for the comparison of TBI patients with their matched controls within homogeneous subgroups.We demonstrated that within each phenotype, patients with TBI are more likely to have NPD than matched controls, corresponding to meaningful differences in NPD risk when matched controls are considered.It is important to note that latent classes were derived using a broader set of indicators, including demographics and injury-related variables for patients with TBI.The matched controls, however, were matched to patients with TBI patients based on age, sex, comorbidities, and psychiatric conditions.This difference should be taken into account when interpreting the results of the conditional logistic regression with each phenotype.
Guidelines from the American College of Surgeons emphasize the importance of postinjury mental health disorder screening and intervention for trauma patients [36].Triaging the risk of NPDs at the time of discharge from hospital may lead to better post-injury mental health outcomes and quality of life for patients with TBI.Our approach of identifying interpretable phenotypes with different risk of NPDs, and the accurate classification of individuals into these phenotypes, allows for a practical risk stratification approach of patients with TBI.The parameter estimates from the LCA model on the original data can be applied to calculate the posterior class membership probabilities for new patients.This process facilitates risk assessment by classifying new patients into one of the defined phenotypes.

Strengths and limitations
This study has several strengths.The linked data sources provided detailed information for a population-based cohort, limiting selection and recall biases, while also ensuring robust statistical power.The variables used in our modelling approach are available in electronic health records, and therefore the insights from our analyses are applicable within similar settings.To assess the health burden of NPDs among individuals with and without TBI, a direct comparison of data between TBI cases and uninjured controls is essential.Hence, we employed a demographically similar uninjured cohort as the matched group for patients with TBI.This matching approach enables the exploration of differences in NPD prevalence between injured individuals and healthy controls.
There are limitations to be considered when interpreting the current findings.Our results are dependent on patient encounters documented in administrative databases, with the potential for under-reporting of diagnoses, lack of specificity in coding, and inaccuracies in designating diagnostic categories.Given that the current study used administrative health databases, we did not capture lifetime psychiatric history using structured clinical interviews.Our study may have potentially underestimated the true prevalence of pre-injury psychiatric history and its impact on the development of novel psychiatric disorders post-TBI using a lookback period of two years to ascertain pre-injury psychiatric burden.For patients with pre-existing psychiatric disorders, we were unable to determine whether these disorders worsened in severity based on ICD coding.Additionally, our study did not capture changes in life satisfaction, functional independence, and social participation, all of which are known to influence mental health functioning after injury [37,38].It is also important to acknowledge that while the control group was selected to match the TBI cohort demographically and did not have a recorded history of TBI or develop TBI during the two years of follow-up, some controls may have experienced head injuries which did not result in clinical encounters.Finally, the current study did not perform external validation of the identified phenotypes.To confirm the generalizability of these phenotypes, future research should apply the posterior membership probabilities derived from our main analysis to a new, independent sample of TBI patients.Nevertheless, the study findings provide a method to identify homogeneous subgroups of TBI patients in relation to NPDs, thereby facilitating both research and clinical applications.We demonstrate that LCA can distinguish clinically meaningful phenotypes using routinely collected variables.Using these standard variables may be a useful method to stratify patients into more homogeneous groups for enrollment into clinical trials, selection of pharmacological and non-pharmacological interventions, and prediction of clinical outcomes.From a clinical standpoint, the risk phenotypes would be useful for counselling patients and predicting their psychiatric rehabilitation needs.

Conclusions
This study highlights the significant post-injury psychiatric burden experienced by survivors of TBI.Head-injured patients with and without a pre-injury psychiatric history have a significantly greater burden of NPDs compared with population-matched controls.TBI is a heterogenous disorder comprised of clinically distinguishable phenotypes with different patterns of NPD risk.Using this approach for risk stratification at the time of discharge from acute care may help guide early and targeted interventions for those individuals most at risk of developing NPDs following TBI.
, motor vehicle collisions (MVC), violence, or other NSTR Injury severity Defined by the maximum Abbreviated Injury Score (AIS).AIS Head between 2-5 were considered as categorical variables NSTR Length of stay Numbers of days at the acute care facility.Categorized as prolonged if LOS ≥ 21 days.NSTR Discharge destinationThe destination of the patient upon discharge from the acute care facility.Categorized as another acute care facility, chronic care facility, nursing home, rehabilitation facility, home with support services, or home.NSTR 1 NSTR: Nova Scotia Trauma Registry 2 HDNS: Health Data Nova Scotia3 Only psychiatric conditions with a prevalence of > 10% were included in LCA.The full list of psychiatric disorder categories is presented in Supplementary Table1

Fig. 1
Fig. 1 Study design.Propensity score matching of TBI vs. control groups (left panel).LCA for two TBI cohorts: those with and without pre-injury psychiatric conditions (right panel)

Fig. 2
Fig. 2 Radar plots for four phenotypes of TBI patients without psychiatric conditions

Table 2
Characteristics of the matched cohort 1All patients within the outcome of NPD-A, did not have a documented history of psychiatric disorders in the two years prior to injury2All patients within the outcome of NPD-P had a documented history of psychiatric disorders in the two years prior to injury3A small proportion of injuries were due other mechanisms including burns, drowning, or asphyxia

Table 3
Results of the conditional logistic regression with adjusted odds ratio (OR) and 95%confidence interval (CI) for NPDs in the matched cohort

Table 4
Fit statistics for latent class analysis The SABIC and LRT indicated the 5-class model as the most appropriate.Lower values for BIC and SABIC indicate relatively better balance between parsimony and model fit.Emphasis was placed on BIC, SABIC, and VLMR-LRT given evidence showing their unique strength in identifying the ideal number of classes.The 4-class model provided a more clearly defined representation of the data over the 5-class models.For these reasons, the 4-class model was selected 5Vuong-Lo-Mendell-Rubin adjusted likelihood ratio test *

Table 5
Analysis of patient characteristics by latent class for the TBI cohort without pre-injury psychiatric conditions.Adjusted relative risk (RR) of NPD-A is estimated from log binomial regression model.Odds ratio (OR) comparing patients with TBI to their matched controls is estimated using conditional regression for matched subjects